Self-sustained emission in semi-infinite non-Hermitian systems at the exceptional 

point 
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Complex potential and non-Hermitian hopping amplitude are building blocks of a non-Hermitian 
quantum network. Appropriate configuration, such as PT-symmetric distribution, can lead to the 
full real spectrum. To investigate the underlying mechanism of this phenomenon, we study the 
phase diagram of a semi-infinite non-Hermitian system. It consists of a finite non-Hermitian cluster 
and a semi-infinite lead. Based on the analysis of the solution of the concrete systems, it is shown 
that it can have the full real spectrum without any requirements on the symmetry and the wave 
function within the lead becomes a unidirectional plane wave at the exceptional point. This universal 
dynamical behavior is demonstrated as the persistent emission and reflectionless absorption of wave 
packets in the typical non-Hermitian systems containing the complex on-site potential and non- 
Hermitian hopping amplitude. 

PACS numbers: 03.65.-w, 11.30.Er, 71.10.Fd 



I. INTRODUCTION 

Non-Hermitian Hamiltonians are often employed to de- 
scribe the open systems due to their features of complex- 
valued energy and non-preserved particle probability. 
Recent observations show that a large families of non- 
Hermitian Hamiltonians can have all eigenvalues real, if 
the loss and gain are set in a balanced manner, being in- 
variant under the combination of the parity (V) and the 
time-reversal (T) symmetry. A parity-time {VT) sym- 
metric non-Hermitian quantum theory has been well de- 
veloped as the complex extension of conventional quan- 
tum mechanics 1-8]. Although the condition of the VT 
symmetry for the complete real spectrum is weaker Q, 
it still implies the underlying mechanism can be based 
on the balance of the loss and gain. However, such an 
intuitive consideration of the balance needs to be investi- 
gated precisely. The concept of the balance should not be 
simply understood as the conjugate relation of two non- 
Hermitian subsystems arising from the VT symmetry. 
It can not provide physical explanation to the following 
features about exceptional point: (i) The VT symmetry 
of the system can not guarantee the balance of the loss 
and the gain, or the reality of the energy levels, (ii) The 
spontaneous symmetry broken states always appear in 
pair. Furthermore, this consideration is also related to 
the precise physical significance of the complex potential 
and non-Hermitian coupling, which are basic elements for 
a discrete non-Hermitian system. On the other hand, the 
purpose of this investigation is not only for the fundamen- 
tal physics, but also for the application in practice due to 
the formal equivalence between the quantum Schrodinger 
equation and the optical wave equation Further- 
more, the VT symmetry breaking has been observed in 
experiments [2ll |22|. 

In this paper we try to clarify the concept of balance 



in the non-Hermitian discrete system in the framework 
of the quantum mechanics rather than the phenomeno- 
logical description. We study the exceptional point of 
a semi-infinite non-Hermitian system from the dynami- 
cal point of view. We will show that the wave function 
within the lead becomes a unidirectional plane wave at 
the exceptional point. This universal dynamical behavior 
is demonstrated as the self-sustained emission and reflec- 
tionless absorption of wave packets by two typical non- 
Hermitian clusters containing the complex on-site poten- 
tial and non-Hermitian hopping amplitude. 

This paper is organized as follows. In Section[H]we an- 
alyze the classification of possible solutions. Section IIIII 
consists of two exactly solvable examples to illustrate our 
main idea. Section ITVl presents the connection to the VT 
symmetric systems. Section [V] is devoted to the numer- 
ical simulation of the wave packet dynamics to demon- 
strate the phenomena of the persistent emission and re- 
flectionless absorption. Section I VII is the summary and 
discussion. 



II. SEMI-INFINITE SYSTEM 



The non-Hermitian tight-binding model, with the non- 
Hermiticity arising from the on-site complex potentials as 
well as the non-Hermitian hopping amplitude, is a nice 
testing ground to study the basic features of the non- 
Hermitian system. In recent years [23l - [28| , it has received 
intensively study not only because of its analytical and 
numerical tractability but also the experimental accessi- 
bility. 
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Let us now consider a non-Hermitian finite cluster cou- 
pled to a semi-infinite lead. The Hamiltonian is written 



2 



as 

H = H t + H sub , 
o 

H l= -jJ2 (afoi+1+H.c), (1) 

/ — — oo 

H sub = ^ ' K ij a i a j- 
ij = l 

It is noted that -ff S ub is non-Hcrmitian, possessing the 
complex- valued eigen energy, while Hi is Hermitian, hav- 
ing complete spectrum E — — 2Jcos(fc), (fee [0, 27r)) 
and the eigen wave function sin (kj) a] |0). To investigate 
the role of the lead in the non- Hermitian H sub , we will 
consider the whole solution of the Hamiltonian H and 
analyze its properties. 

The wave function can be expressed as \ip) = 
J2j=-oo f k (J) a ] |0)- The explicit form of the wave func- 
tion f k (j) depends on the structure of -ff su b- General 
speaking, the solution of f k (j) can not be obtained ex- 
actly even the explicit form of i? su b is given. However, 
within the lead the wave function is always in the form 

f k (j <0) = A k e ik * + B k e~ ik ^ (2) 




(b) 

FIG. 1. (Color online) (a) Schematic illustration of the con- 
figuration of the concerned network and profiles of possible 
solutions. It consists of a non-Hermitian finite cluster and a 
semi-infinite chain. The solutions on the lead are threes types: 
scattering state wave (green), monotonic damping wave (blue) 
and oscillation damping wave (red), (b) When an arbitrary 
semi-infinite system and its conjugate counterpart are con- 
nected at infinity, a "PT-symmetric system are constructed. 



due to the semi-infinite boundary condition. 
The Schrodinger equation has the explicit form 

- Jf k U - 1) - Jf k (J + 1) = Ef k (j) , (j < 0) 

-Jf k (0) + E^^iif k ii)=Ef k a), (3) 

Ef=i^f k H)=Ef k (j), C?e[2,iV s ]) 

within all the regions. The solutions of E and f k (j) de- 
pend on the structure of the system -ff S ub- Nevertheless, 
the exclusive geometry of the lead will give some clues to 
the characteristics of the eigenvalues and eigenfunctions. 

In the framework of Bethe ansatz method, all possi- 
ble solutions within the lead can be classified into three 
types: 

(1) Scattering state wave, the wave function and energy 
are in the form 

/ss (j) = A k e ikj + B k e~ ikj , (k is real) (4) 
E = -2Jcos(fc). (5) 

(2) Monotonic damping wave, the wave function and 
energy are in the form 

/ MD (j) = (±l) j e« CS<0) (6) 
E = -2Jcosh(/3) . (7) 

(3) Oscillation damping wave, the wave function and 
energy are in the form 

/od (j) = e^ + « , 08 < 0, k € (0, 2vr) , k ? tt) (8) 
E = -2 J [cos (k) cosh (/?) + i sin (k) sinh (/3)] .(9) 



In Fig. Q] (a) , the concerned system and three types of 
possible solutions within the lead are illustrated schemat- 
ically. For the case of hermitian ff su b , the solution are the 
form of fss (j) and /md (j) with \A\ or \B\ — definitely. 
In the case of non-Hermitian i? su b, /od (j) may appear 
associated with the complex energy level. In case of ab- 
sence of the solution /od (j), full real spectrum achieves, 
which shows the existence of the stationary states. It 
indicates that the lead acts as a channel to balance the 
gain or loss in the system H sub - 

At certain points k Cl the system makes transitions be- 
tween eigenstates fss (j) and /md (j), as well as between 
/ss (J) and /od (j)- The former transition is actually a 
switch between real and imaginary k, preserving the real- 
ity of the eigen energy. Then the transition point locates 
at k = 0, tt, i.e., 

/md U) ^ (±iy , (io) 

ZssO) —4 (±iy, (ii) 

which usually occurs in the case of Hermitian H sub . The 
later transition only occurs in a non-Hermitian system, 
eigen energy switching between real and complex values. 
In contrast to the above case, the transition point (re- 
ferred as exceptional point) depends on the structure of 
the non-Hermitian H sub , i.e., 

/odCO^H (12) 
/ss (?) ^ e ifc < (13) 

It indicates that a unidirectional plane wave exists in the 
lead when an appropriate non-Hermitian _ff su b is con- 



nected. It has both fundamental as well as practical im- 
plications. This result reveals the exceptional point from 
an alternative way: It is the threshold of the balance be- 
tween the non-Hcrmitian subcluster and the lead. From 
a practical perspective, the unidirectional-plane- wave so- 
lution at the exceptional point can be used to realize 
the reflectionless absorption and persistent emission in 
the experiment. To characterize the probability genera- 
tion (negative in the case of the dissipation) of the non- 
Hermitian cluster, we introduce the current operator 



-iJ 



c j+l c j 



(14) 



where j e [—00, 0]. 

For three types of eigenstates fss(j) e ~ lEt : 
/md (j) e~ lEt and /od (j) e~ tEt , the corresponding 
currents can be obtained as 



J ss = 2J(\A k \ 2 -\B k \ 2 )sm(k), 
Jmd = 0, 

Job = 2Jsin(fc) e 2 [«- 2Jsin W sin ^'] 



(15) 



We can see that Jss is time-independent and is conser- 
vative along the lead, representing a steady flow or the 
dynamic balance, while Jod is non-periodically time- 
dependent, indicating the unbalance of the state. In 
other word, the mechanism of the reality of the spectrum 
is the balance between the source (or drain) and the chan- 
nel of the probability flow. Then the exceptional point is 
the threshold of such dynamic balance, corresponding to 
the unidirectional-plane- wave, i.e., \A k \ — or \B k \ = 0. 
Then the probability generation for the exceptional point 
is 



J c = 2Jsin (fc c ) 



(16) 



which the sign indicates that the cluster is a source or 
drain, then it is referred as critical current in this paper. 
Unlike the situation in traditional quantum mechanics, 
the magnitude of the current J c does not represent the 
absolute current in traditional quantum mechanics be- 
cause the corresponding eigenstate is not normalized un- 
der the Dirac inner product. 

We will demonstrate and explain these points through 
the following illustrative example. We would like to point 
out that, there is an another type of the exceptional 
point, arising from the transition of two types of bound 
states /od (j) and /md (j), which is beyond our interest. 



III. ILLUSTRATIVE EXAMPLES 

In this section, we investigate two simple exactly solv- 
able systems to illustrate the main idea of this paper. In 
order to exemplify the above mentioned analysis of relat- 
ing the wavefunction within the lead and the eigenvalue, 
we take iJ su b to be the simplest non-Hermitian networks 
to construct two types of exemplified systems. Type I is 
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FIG. 2. (Color online) Schematic illustration for the non- 
Hermitian 'PT-symmetric networks: (a) an uniform chain 
with a complex potential and (b) with a complex hopping 
at both sides. 



a uniform chain with a complex potential at one end and 
type II is a uniform chain with a complex hopping at one 
end. In Fig. [5J the ending complex potential and hop- 
ping are illustrated schematically. In the following, we 
present the analytical results in the framework of above 
mentioned for the two models in order to perform a com- 
prehensive study. 



A. Complex potential 

The type I Hamiltonian has the form 

00 

Hcp = -J^(a]a J+ i + H.c.) + Je ie ala , (17) 
3=0 

where 9 is a complex number. According to Bethe ansatz 
method, the wavefunction f k (j) can be expressed as 

f k (3) = A k e^ + B k e- ik \j € [0,oo) (18) 
and the Schrodinger equations for Hqp is 

- J.f k (3 + 1) - Jf" (3 - 1) = E k f k (j) , (19) 
-Jf k (1) = [E k - Je w ] f k (0) . 

Submitting f k (j) into the Schrodinger equation, we have 



D K 



B k e~ lk + c 



(20) 



which is the reflection amplitude for the scattering state. 
Now we are interested in the wavefunction with complex 
eigen energy. The existence of the solution /od (j) re- 
quires 



-e l8 and Im (fc) > 0, 



(21) 
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FIG. 3. (Color online) Schematic illustration of the phase 
diagrams in the infinite N system: (a) a uniform chain with 
a complex potential and (b) with a complex hopping at one 
end. 



which leads k — it — 9 with lm(9) < 0. Then we conclude 
that there is a unique complex solution within the region 
lm(0) < and the system has full real spectrum if the 
potential is in the rest region. At the boundary, we have 



Re (0 C ) = 6 C 



(22) 



which indicates a circle of radius J in the complex plane. 
The phase diagram is sketched in Fig. [3] (a). Then the 
corresponding wavefunction has the form 

f K (j) = e«*- e <>», (23) 
which represents a unidirectional plane wave with energy 

E kc = 2Jcos(6» c ) . (24) 
Accordingly, the critical current is 

J c = 2Jsm(6 c ), (25) 

which accords with the intuition that a positive imagi- 
nary potential can be a source and a negative imaginary 
potential can be a drain. However, unlike the situation 
in traditional quantum mechanics, the magnitude of the 
current J c does not represent the absolute current in tra- 
ditional quantum mechanics because the eigenstate is not 
normalized under the Dirac inner product. Before further 
discussion of the implication of the obtained result, two 
distinguishing features need to be mentioned. Firstly, the 



non-Hermitian system can have full real spectrum even 
though there is no symmetry required. Secondly, there 
is only one possible complex energy level. Then there 
is no level coalescing occurring at the exceptional point. 
Both of these differ from that of a finite 'PT-symmetric 
system. In the following example, it will be shown that 
such features are not exclusive to the complex potential. 



B. Complex-coupling dimer 

The type II Hamiltonian has the form 

oo 



'CC 



j\^(atoj+i +H.c.) — Jna^ai — Jn*a\aa, (26) 

3=0 



which is a different type of the non-Hermitian model in 
contrast to the type I. Here k is a complex number. In the 
following, we will perform a parallel investigation with 
the current Hamiltonian. The Bethe ansatz wave function 
has the form 



f k U) 



A k e^+B k e~^, j G [l,oo) 



C k e ikj + D k e~ ik ^ j = 
Substituting f k (j) to the Schrddinger equation, 

- Jf U - 1) - Jf U + 1) = E k f k (j) , j G [2, oo) 



(27) 



-J K f k (0) - Jf k (2) = E k .f k (1) , 
-Jnf k (1) = E k f k (0) , 

we obtain the reflection amplitude 

B k (k 2 - 1) e lk 



- ik 



R k 



A k (k 2 - 1) e~ ik - e ik ' 
The existence of the solution /od (J) requires 

k = -iln (k 2 - l) and Im (k) > 0. 



(28) 



(29) 



(30) 



Similarly, we conclude that there is a unique complex 
solution within the region \k | - 2Re(K 2 ) > 0, and the 
system has full real spectrum if the potential is in the 
region. The phase diagram is sketched in Fig. [3] (b). 
The boundary can be expressed as 

i { [Re (k c )} 2 + [Im (k c )] 2 } 2 = [Re ( Kc )] 2 - [Im (k c )] 2 , 

(31) 

which is a Lemniscate of Bernoulli equation in the com- 
plex plane. Then the corresponding wave function has 
the form 

f*-m-/ e ^ /2b ' 3 e [i,oo) ( 

f U) - i K-WVM, j = ' (32) 



where 



tan ip 



(33) 



5 



is real. The wave function f ka (j) represents a unidirec- 
tional plane wave with energy Ek c = — 2Jcos(<^/2) and 
the critical current is J c = 2Jsin(ip/2). It can be shown 
by straightforward algebra that J c > (,7c < 0) in the 
second and fourth (the first and third) quadrants in the 
complex K-plane. It indicates that, like the complex po- 
tential, a complex-coupling dimer can also be a source or 
drain. 

It is noted that the systems in both above two exam- 
ples are not symmetric, which show that the symmetry 
is not the necessary condition for the occurrence of full 
real spectrum. The underlying mechanism can be ex- 
plained as the balance between the source (or drain) and 
the channel. In this sense, a semi-infinite chain can act 
as a source (or drain) to balance the original drain (or 
source). The single exceptional point is the threshold of 
such balance. This point will be elucidated in detail in 
the section HVl 



IV. CONNECTION TO VT SYMMETRIC 
SYSTEMS 

Now we consider the connection between the ob- 
tained results and the well developed non-Hermitian 
■pT-symmetric quantum mechanics. Intuitively, a VT- 
symmetric system can be simply constructed by connect- 
ing an arbitrary semi-infinite system and its conjugate 
counterpart at infinity, as sketched in Fig. [Tj (b). We 
will demonstrate this point through the following illus- 
trative examples. We start from the corresponding solu- 
tions of the Hamiltonians Hq P and Hq C , which actually 
can be obtained by applying time-reversal operation, i.e., 
taking complex conjugation for the obtained solutions of 
the Hamiltonians i?cp an d Hqc- It is interesting to find 
that the real eigen-valued solutions between Hqp and 
H^p (Hqc and H^ c ) can match with each other due to 
the fact that 



Ak 
B k 



A _» 



(34) 



for both examples in the Eqs. (|2T)]) and (121)1) . In other 
word, all the real eigen-valued solutions will not change 
if two systems £/cp and H^ p (Hqc and H^ c ) are con- 
nected at the infinity. Although this is not a surprising 
result, we still verify it explicitly in a strict manner by 
solving the VT non-Hermitian systems, which are finite 
versions of combined -ffcp and Hq P (Hqc and Hq C ). A 
sketch of such systems are given in Figs. [5] (a) and [5] (b). 

The first example is a VT symmetric non-Hermitian 
TV-site chain with complex on-site potential at two ends, 
which has the Hamiltonian 

N-l 

H% = — J (atctj+i + H.c.) + Je ie a\ax + Je~ ie * a* N a N - 

(35) 

It is a VT symmetric model, i.e., [VT,Hi] = 0, where 
the action of the parity operator V is defined as V : I — » 



N + 1 — 1 and the time-reversal operator T as T : i — > — i. 
For infinite N, it becomes the combination of the systems 
Hqp and Hq P . For finite N, it is an extension version of 
the model proposed in the previous paper Ref. (23j. By 
using the standard Bcthc ansatz method, the solution is 
determined by the critical equation 



r(jfc) =0 



(36) 



where 



r(fc) 



" > sin [k(N - 1)] (37) 
- sin [k(N + 1)] + 2Re (e l9 ) sin (kN) . 



Accordingly, the exceptional point can be obtained by 
the equations [H, 



T (fc c ) = and 



dT(fc) 



dk 



0. 



(38) 



k=k c 



It is difficult to get the explicit solutions of the Eq. 
for finite N. Nevertheless, the equation about dr (k) /dk 
can be reduced to 

e i(e " r) cos [k c (N - 1)] + cos [k c (N + 1)] (39) 
+2Rc (e ie ) cos (k c N) 

by taking the approximation N ± 1 N in the large 
N limit. It is easy to find that the existence of real k c 
solution requires Im((9) = and k c = tt — 9. It is in 
accordance with result in the corresponding semi-infinite 
system. 

The second example can be described by the Hamilto- 
nian 



N-2 



H 



— J (ajaj-|_i+H.c.)— J(Ka[a2+K*a J N _ 1 aN+H..c.), 

(40) 

which corresponds to the combination of two systems 
Hqc and H^ c . By the same procedure as that for Hi, 
we find that the critical equation for finite N is 



X2sm[k (N — 3)] + \i sin [k (N — 1)] 
- sin [k (N + 1)] = 



where 



Xi 

X2 



K 2 + (K*Y 

xi - M 4 - 



l. 



(41) 



(42) 
(43) 



In addition, the exceptional point for large N is deter- 
mined by the equations 



and 



2 X 2 cos (2fc c ) + xi « 



2cos(2fc c ) -xi = 0, 



(44) 



(45) 



which leads to the same results as the Eqs. (|3~Tj) and 
The above results are helpful to understand the 
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FIG. 4. (Color online)The profile of the evolved wave function {j \4>(t)) for the initial state being (a) an incident Gaussian 
wavepacket with ko = — n/2, Na = 10 and a — 0.5, and (b) a site state at the left end. We plot the imaginary (blue) and real 
(red) parts of the amplitudes {j \<j> (t)) of the wave function at instants t with the unit of 1/J. We see that in both cases, the 
evolved wave functions tend to the persistent emission of the plane wave with the momentum k ~ k c — n/4. One can see the 
plane wave travels to the left and the corresponding phase velocity is 4\/2/tv from the subfigures t — 80 and t = 80.5 (t = 60 
and t = 60.5), in the Figs. g](a) 01(b)). 



mechanism of the Hermiticity of a non-Hcrmitian system. 
The existence of the full real spectra of two above VT 
systems is attributed to the balance between the source 
and drain. The accordance between the results of these 
VT systems and that of semi-infinite systems implies that 
a lead can act as a multifunctional source (or drain) to 
balance any arbitrary drain (or source). 



V. WAVEPACKET DYNAMICS 

In this section we will apply the above theoretical re- 
sults to simple accessible examples to investigate the dy- 
namic behavior for local initial states. This may provide 
some insights into the application in practice. 

Firstly, we focus on the phenomenon of the persistent 
emission. Consider an arbitrary local initial state on the 
lead in the system at the exceptional point. The initial 
state can always be written in the form 

(j \</> (t = 0)) = c kc e ik ^ + $> fc ( e ' ifcJ + R k^ tkl ) , 

( 46 ) 

where ^2 k ^ k c k (e 1 *^ + R k e represents the superpo- 
sition of the scattering states with different k. It is pre- 
sumable that the probability of all the scattering states 
transfers to infinity after a sufficient long time, and then 
only the unidirectional plane wave survives. To demon- 
strate and verify this analysis, numerical simulations are 
performed for two typical initial states: an incoming 
wavepacket and a site state at the scattering center. A 
Gaussian wavepacket with momentum k and initial cen- 
ter Na has the form 

\i>(k ,N A )) = J=Y,e-^ U - NA)2 e lk ° 3 \J) , (47) 



where flo = 52 ■ e a ^ Na ^ is the normalization factor 
and the half- width of the wavepacket is 2 /a. Here 
we take fcrj = ~ ' r /2, Na = 100 and a = 0.5. The con- 
cerned system is described by the Hamiltonian in the Eq. 
(|35|) with 9 — 37r/4, which corresponds to the persistent 
emission of the plane wave with momentum k c = 7r/4. 
The profiles of the evolved wave functions are plotted in 
Figs. |4] (a) and|4] (b). One can see the profile of the wave 
and the corresponding phase velocity from the figures, 
and after a little long time the evolved wave functions 
accord with the plane wave of 

Re ((j |0 (*))) ~ cos [(tt/4) j + 2Jt cos (tt/4)] , (48) 
Im((j |0(t)))~ S m[(7r/4)i + 2Jtcos(7r/4)], (49) 

within the finite region along the lead, approximately. 
This result has implications in two aspects: First, we 
achieve a better understanding of the imaginary poten- 
tial. We found that a complex potential always corre- 
sponds to the wave vector of the unidirectional plane 
wave, which is determined by the Eq. (|2"5)l . Second, 
it provides a method to measure the complex potential 
in the experiment. 

On the other hand, it is presumable that the reflec- 
tionless absorbtion should naturally be reflected in the 
dynamics of the wavepacket with the momenta around 
k c due to the continuity of the reflection coefficient in 
the vicinity of the exceptional point k c . Consider an in- 
coming Gaussian wavepacket with momentum ko = k c 
and initial center Na, which can always be written as 



(j \i>(k c> N A )) = 




(50) 

where (j | k) = e %k i denotes the plane wave with momen- 
tum k, and = f-e~( k ~ kc ' ^ 2a dfc is the normalization 
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FIG. 5. (Color online)The probability distribution of the evolved wave function for the initial state being the incident Gaussian 
wavepackets with ko = — vr/2, Na = 100, (a) a = 0.5, and (b) a = 0.015. We plot \(j \cf)(t))\ 2 at instants t with the unit of 
1/J. We see that in both cases, the reflection coefficients are small, especially for the wider incident wavepacket, which is more 
close to the plane wave with the momentum k — k c — tt/2. 



factor. It is a superposition of the plane waves with mo- 
menta around k c , which have small reflection coefficients. 
Therefore there is small probability being reflected for the 
wavepacket. The optimal situations to achieve the reflec- 
tionless absorbtion of a wavpacket are determined by the 
conditions 



oretical analysis. 



i? fe J = 0, 



d\R k \ 



dk 



= 0, 



(51) 



which minimizes the reflectional probability. Straight- 
forward derivation indicates that the optimal complex- 
potential (hopping) scattering center system requires 9 = 
tt/2 (k 2 = 1 ±i) and the momentum of corresponding re- 
flectionless wave is k c = tt/2 (tt/4). 

To demonstrate and verify this analysis, numerical sim- 
ulations are performed for the two initial wavepackets 
with k = -tt/2, N a = 100, a = 0.5 and 0.15, re- 
spectively. The concerned system is described by the 
Hamiltonian in the Eq. (|35|) with 9 = tt/2, which corre- 
sponds to the reflectionless plane wave with momentum 
k c = tt/2. The profiles of the evolved wave functions are 
plotted in Figs. [5] (a) and (b). The reflection coefficients 
of the two wavepackets are 0.003 and 0.035, respectively. 
It shows that the wider of the wavepacket the lower of 
the reflection rate, which accords with the previous the- 



VI. SUMMARY 



In summary, the mechanism of the non-Hermiticity of 
a discrete non-Hermitian system has been investigated 
in an alternative way. It is shown that the symmetry 
is not the necessary condition for the occurrence of full 
real spectrum. The underlying mechanism can be ex- 
plained as the balance between a non-Hermitian cluster 
and a semi-infinite chain: the semi-infinite lead can play a 
complete role to balance a finite non-Hermitian cluster, 
resulting full real spectrum and the single exceptional 
point as the threshold of the such balance. Furthermore, 
at the exceptional point, the eigen wave function is shown 
to be a unidirectional plane wave. Practical application 
of this feature to the dynamics of the wave packet demon- 
strates the phenomena of the self-sustained emission and 
reflectionless absorbtion, which could be applied to the 
quantum device design. 
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